Cooper-pair resonances and subgap Coulomb blockade in a superconducting 

single-electron transistor 
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We have fabricated and measured superconducting single-electron transistors with Al leads and 
Nb islands. At bias voltages below the gap of Nb we observe clear signatures of resonant tunneling 
of Cooper pairs, and of Coulomb blockade of the subgap currents due to linewidth broadening of 
the energy levels in the superconducting density of states of Nb. The experimental results are in 
good agreement with numerical simulations. 
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The single-electron transistor— (SET) and its supercon- 
ducting version is one of the most versatile tools in meso- 
scopic physics. It has been used for extremely sensitive 
charge measurements)^ for the construction of Cooper 
pair pumps and other adiabatic devices with appHcations 
in metrology, and more recently for building up super- 
conducting quantum bits 4 

The IV characteristics of superconducting SETs 
present the usual features of quasiparticle tunneling (at 
voltages above 2ANb + 2Aai), Josephson-quasiparticle 
tunneling (at half of these values), and Josephson effect 
(around zero bias). These features have been thoroughly 
investigated by now by many groups and the physics of 
a charge transport at these bias voltages is well under- 
stood. However, at low bias voltages also other trans- 
port processes could become important and can alter the 
performance of Josephson-based devices. In this paper, 
we study two such processes appearing in our Nb-based 
SET: Resonant tunneling of Cooper pairs, and transport 
through states inside the gap of Nb (subgap currents). 

We have fabricated Al/AlOa^/Nb/AlOa^/Al single elec- 
tron transistors using a Uthographic technique described 
elsewhere Measurements were done using a small dilu- 
tion refrigerator equipped with well-thermalized and elec- 
trically filtered measuring lines. The superconducting 
gaps obtained for Nb and Al (Awb = 1.4 mV, Aai — 0.2 
meV), and also the measured critical temperatures for 
Nb (Tc,Nb ~ 8.0 — 8.5 K), show that the films are indeed 
of good quality. 

At voltages below the gap of Nb, a series of gate- 
dependent resonance peaks appears in the IV charac- 
teristics of the SET (Fig. [T]). We interpret this as reso- 
nant tunneling of Cooper pairs, a transport phenomenon 
first predicted theoretically and later observed in Al 
symmetrically-biased superconducting SETsj^ Below we 
describe the same process for our Nb-island SETs un- 
der the asymmetric bias shown in Fig. [TJ We consider a 
generic process in which a charge 6qi tunnels through the 
left junction and a charge Sq2 tunnels through the sec- 
ond junction, both into the island. During the process, 
a charge Sq =^ 5qi + Sq2 is transferred into the island 
and a charge 5Q ~ 5qi — 5q2 is transferred through the 



external circuit in the forward direction. The change in 
the electrostatic free energy (including work done by the 
sources) associated with this process is 
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where go is the initial charge of the island, Ci and C2 
are the capacitances of the left and right junctions, Cg is 
the gate capacitance, and Cs = Ci -|- C2 + Cg . Resonant 
Cooper pair tunneling in superconducting SETs occurs 
when no energy is required for processes resulting in the 
transport of m Cooper pairs SQ = —2'me through the 
external circuit and the creation of an excess of n Cooper 
pairs Sq = — 2rte on the island, i.e., E{—2ne, — 2me) = 0. 

The dominant processes are those involving m — ±1 
and n = ±1, corresponding to a single Cooper pair tun- 
neling through either one of the junctionsi^ For these 
processes, the resonant condition describes two lattices 
with a cell size of (4e/Cs)x(2e/Cg) (in VxVg plane) cor- 
responding to odd and even values of qo/e, and displaced 
with respect to each other by 2e /Ce along the V axis and 
by e/Cg along the the Vg axis. During the time of an ex- 
periment, qo can either be fixed, in which case only one 
lattice pattern will appear, or it can fluctuate between 
odd and even values, in which case what will be mea- 
sured is the overlap between the two lattices, resulting in 
a checkerboard lattice with half the periodicity. This last 
situation occurs indeed in all of our three samples. To 
check this, we first measured the gate modulation at volt- 
ages above the quasiparticle threshold of Nb, where the 
transport is of single-electron type, and determined e/Cg. 
This corresponds to the size of the checkerboard pattern 
along Vg that we see in our measurements at lower bias 
as wen, e.g., e/Cg = 12.3 mV (Cg = 12.9 aF) for the 
sample of Fig. [1] (For clarity, we will concentrate our 
discussion on that particular sample from now on. The 
other two samples yielded similar results.) The observa- 
tion that the checkerboard pattern is le-periodic allows 
us to determine e/Cs = 68.3 /xV and thus the charging 
energy Ec = eV2Cs = 34.1 i^eVJ- The SET's junction 
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asymmetry is reflected in the difi^erent absolute values of signs) versus the negative-slope lines (m and n have the 
the positive-slope lines (for which n and m have different same sign), yielding Ci = 1.03 fF and C2 — 1.29 fF. 




FIG. 1: Schematic of the ch-cuit (upper left inset) and current versus bias and gate voltages (3D plot), showing the checkerboard 
pattern (right inset) associated with resonant tunneling of Cooper pairs. The contour plot of the right inset is obtained by 
subtracting the background current thus leaving only the characteristics of the resonant tunneling. 



To check for consistency, we have determined Eq also 
by standard Coulomb-blockade thermometry^ with the 
superconductivity both in Al and Nb suppressed. We 
find Eq ~ 36 /ieV, which is in very close agreement with 
the results above. 

The second intriguing feature in the low-bias IV of our 
Nb-island SETs, is the presence of relatively large sub- 
gap currents: from figures[T]and[Ha) one can see that the 
Cooper pair resonances appear as being built upon a cur- 
rent which increases with bias voltage (at voltages larger 
than 2Aai -I- 2ANb, it will merge into the quasiparticle 
current). The origin of the subgap currents in supercon- 
ductors is still a matter of intense theoretical and exper- 
imental investigations.— Here, we further explore these 
currents by using a magnetic field to suppress the gap of 
Al; the transition to the normal-metal state in the leads 
is clearly indicated by the disappearance of the Joseph- 
son effect and of the Cooper pair resonances as seen in 
Fig. Ha). 

While Al is in the normal state, a very interesting fea- 
ture appears at low bias voltages [Fig. [2Kb)] , where we ob- 
serve a dip in the conductance. In measurements on sin- 
gle Al-Nb junctions [again with the gap of Al suppressed, 



Fig. [2Kb) inset] this feature is not present Instead, below 
0.2 mV the IV is approximately ohmic. This shows that 
the dip is due to the Coulomb blockade of the subgap 
current, and that in the first approximation one could 
attempt to fit the graph with the standard analytical 
expressions for Coulomb-blockade thermometry* (CBT) 
[Fig. [2Kb), main graph, red dashed line]. Fitting yields 
T « 98 mK, which agrees well with the value measured 
using a calibrated resistor thermally anchored to the mix- 
ing chamber, and Eq ^ 13 /xeV, which is in reason- 
able agreement with the value determined from Cooper 
pair resonances, considering that we are not in the limit 
k^T ^ Ec and therefore the standard approximations 
of Coulomb blockade thermometry are not accurate, es- 
pecially for determining Eq-' 



A better model should take into account that we are at 
low temperatures k^T <Ec and also describe the nonlin- 
ear increase in the subgap current at higher voltages. To 
develop such a model, we introduce a lifetime broadening 
r of the quasiparticle energies, resulting in a density of 
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FIG. 2: (a) Effect of a magnetic Geld on IV characteristics 
(current values at different Vg 's, showing gate modulation, are 
superposed): zero, upper curve; intermediate value, middle 
curve; complete suppression of the gap of Al, lower curve, (b) 
The main graph shows the low-bias Coulomb blockade dip: 
experimental values (blue), standard BCS Coulomb-blockade 
thermometry fit (red dotted line), and the predictions of the 
orthodox theory with lifetime broadening and gap inhomo- 
geneity (black continuous line). The upper inset shows the 
conductances for two single Al-Nb junctions of normal resis- 
tances 30 kfl (lower curve) and 11 kQ (upper curve). 
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This effect could be caused by the proximity to the metal- 
insulator transition, due to the granular structure of the 
Nb films (with our films being still on the metalic side, 
r <C ANb)*^ Another possibility is the opening of con- 
duction channels in the junctions^ as suggested by the 
relatively large value of the room-temperature conduc- 
tance [for the sample presented in Fig. [TJ Fig. [21 and 
Fig. 3(a), this was 30 ^S]. Irrespective to the microscopic 
origin, this density of states accounts well for the ex- 
istence of subgap currents seen in Fig. [2l^b), inset. At 
low energies E <C the density of states can be ap- 
proximated as p{E,ANb) ~ r/ANb + (3r/2A)(i;/ANb)^- 



In the case of a single junction of resistance R with a 
normal-state metal this results in a subgap conductance 
i?~^p(y, ANb); therefore, ohmic for low enough voltages 
and increasing as for relatively larger voltages. In 
addition, the granular structure of the films and, pos- 
sibly, impurities resulting from outgassing of the mask 
polymer is also making the film inhomogenous. As a re- 
sult, the gap edge is smeared due to local fluctuations in 
the effective electron-electron interaction. Other fabrica- 
tion techniques [e.g. the use of stronger polymers such 
as PES^^] could improve the quality of the Nb films. For 
our samples, since the superconducting coherence length 
in bulk Nb is only ^Nb = 38 nm, we expect to be in the 
limit in which the size of the inhomogeneities is much 
larger than ^Nb- In this case, it has been shown -'^^ that 
the effective (average) density of states takes the form 
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where V is the probability density associated with a cer- 
tain value of the gap. For small levels of inhomogeneity, 
the function P is a Gaussian of standard distribution a 
centered around an average gap value (slightly smaller 
than the bulk value). 

Finally, to complete our model, we calculate numer- 
ically the currents and conductances for the supercon- 
ducting SET in the framework of the orthodox theory of 
sequential tunneling, by computing the tunneling prob- 
abilities into and out of the island (with the density of 
states above for Nb) and by solving the master equa- 
tion for the charge on the island. In the case of sup- 
pressed Al gap, a very good fit is obtained for a — 0.38 
meV, r — 37.5 ^eV, T — 95.7 mK, and charging en- 
ergy, Ec — 34.1 peV, as determined before [Fig. [U^b), 
main graph, black continuous line]. The temperature also 
agrees well with the value obtained by simple CBT fitting 
above and the value of the thermometer. 
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FIG. 3: (a) Qualitative comparison between the conductances 
extracted from the data by averaging out the gate modulation 
(symbols) and those obtained by numerical simulations (lines) 
at different magnetic helds. (b) The gap of Al leads (open 
circle), measured separately, plotted together with the Al gaps 
determined from the conductance peaks of (a) for three SET 
samples (square, triangle, and star). 

In the case when the leads are superconducting the 
comparison with the experimental data can be done only 
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qualitatively: the main reason is the existence of Cooper 
pair resonances peaks, yielding an extra contribution to 
the current which cannot be eliminated in a straightfor- 
ward way. However, a number of qualitative features can 
still be observed by averaging the resonance peaks over 
gate voltages. The conductances obtained by this proce- 
dure should reflect the voltage dependence of the subgap 
current in the region where the height of the resonance 
peaks is approximately constant. A reasonably fair qual- 
itative agreement with the model presented above is ob- 
tained [Fig. Ela)]. The main feature that we see both 
in the calculated conductivities and in the conductivities 
obtained from our data through the above procedure is 
the appearance of a peak at the onset of the Al quasi- 
particle threshold (2Aai). By using a magnetic field to 
partially suppress the gap of Al, we see that the peak 
moves to the left [see Fig. [2ja) and Fig.IS^a)]. To verify 
that this is indeed the case, we have measured single Al- 
Al junctions with the same fabrication parameters as the 
leads of the Nb- island SET, and determined the gap of Al 
at various magnetic fields [Fig.[3]^b), circles]. The agree- 
ment with the gaps determined from the low-voltage fea- 
tures [Fig.[3Ia)] of three Nb-island SET samples (square. 



triangle, star) is good. 

In conclusion, we have fabricated and measured super- 
conducting single-electron transistors and we have pre- 
sented evidence for a number of transport processes oc- 
curring at bias voltages below the gap of the island: res- 
onant tunneling of Cooper pairs. Coulomb blockade of 
the subgap current for normal leads, and the appearance 
of a step in the subgap current at 2Aai for the case of 
superconducting leads. 
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